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7.  Pgs.  24  and  25,  caption  for  tables: 

Q*1/^1  instead  of  Q^/Q^ 
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10.  P.  33,  equation  following  eq.  (100):  R(x,  x!,P^')  instead  of  R(x,x*  ) 
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VIBRATION  IN  AN  INCOMPRESSIBLE!  FLUID 

Introduction 

In  the  current  procedure  for  calculating  the  natural  frequencies 
of  a  vibrating  ship  the  effect  of  the  fluid  is  taken  into  account  by  increas¬ 
ing  -the  mass  at  each  section  of  the  ship  along  its  length  by  an  added-mase. 

The  latter  is  obtained  by  the  so-called  method  of  strip  theory,  first  pro¬ 
posed  independently  by  F.  M.  Lewis  [1]*  and  J.  L.  Taylor  [2],  in  which  the 
added  mass  at  a  transverse  section  is  first  taken  to  be  that  for  a  two-dimen¬ 
sional  form  of  that  shape,  and  then  the  ordinates  of  this  added-mass  distribu¬ 
tion  curve  are  multiplied  by  a  constant  correction  factor  obtained  from  the 
exact  potential- flow  solution  for  a  vibrating  spheroid. 

The  aforementioned  method  has  yielded  predictions  of  vibration  fre¬ 
quency  in  good  agreement  with  experiment  at  the  lower  modes,  but  serious  devi¬ 
ations  observed  at  the  higher  ones  raise  the  question  whether  more  exact  the¬ 
ories  of  both  the  elastic  characteristics  of  the  ship  and  the  effect  of  the 
surrounding  fluid  might'  not  be  required.  It  lias  been  reported  [5]  that  an 
attempt  is  under  way  to  improve  upon  the  beam  theory  of  ship  vibration  by  tak¬ 
ing  into  account  the  three-dimensional  structure  of  the  hull,  which  clearly 
becomes  important  at  the  higher  modes  of  vibration.  Complementary  to  this  is 
the  work  reported  here  to  improve  upon  strip  theory  by  developing  a  unified 
hydroelastic  theory  of  ship  vibration. 

This  work  was  supported  by  Contract  Nonr  3271(  01 )  { X)  with  the 
Structural  Mechanics  laboratory  of  the  David  Taylor  Model  Basin. 

General  Theory 

First  suppose  that  the  body  consists  of  n  discrete  masses  con» 
neeted  by  certain  elastic  restraints .  Let  q^.  denote  a  set  of  generalized 
coordinates  which  describe  the  displacements,  assumed  small,  of  the  system 
from  equilibrium,  and  q^  the  corresponding  set  of  generalized  velocities. 

'lumbers  in.  [  3  indicate  references  at  end.  of  report. 


Then  the  kinetic  energy  of  the  body,  TB,  is  expressible  in  the  form 
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I,]  =  1,2,  ...  N  „  //=  3n,  (1) 


in  which  the  are  generalized  masses.  In  general  the  are  func¬ 

tions  of  position,  but  for  small  oscillations  about  equilibrium  they  may  be 
assumed  to  have  the  constant  values  corresponding  to  the  equilibrium  posi¬ 
tion,  q,  *0  .  Simultaneously  the  potential  energy  V  of  the  system  may 
be  expressed  in  the  form 


(2) 


in  which,  for  the  same  reason,  the 


may  be  considered  to  be  constants. 


If  the  vibrating  body  is  immersed  in  an  incompressible  fluid,  the 
fluid  will  also  be  set  into  motion  and  have  kinetic  energy  Tf  which  is  ex¬ 
pressible  as  the  quadratic  form  [4] 


2Tf  =  £ 
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L,j  *  1,2,  ...  N 


(3) 


in  which  the  added  masses  a. 


are  also  constants.  The  determination  of 


the  is  a  hydrodynamic,  potential- flow  problem.  This  may  be  approached 

by  introducing  velocity  potentials  0^  corresponding  to  unit  magnitude  of 
each  generalized  velocity  q^  and  formulating  N  Neumann  problems  for  the 
0^  .  As  is  well  known,  the  are  expressible  in  terms  of  the  0^  by 

the  relations  [4] 
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in  which  the  integral  extends  over  the  wetted  surfaces  of  the  system.  By 
formulating  these  Neumann  problems  as  integral  equations  and  replacing  the 
integrals  by  quadrature  formulas,  the  solution  of  the  Neumann  problems  and 
the  evaluation  of  the  are  reduced  to  linear  algebra.  This  procedure 
becomes  feasible  if  high  speed  computing  equipment  is  available. 
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The  dynamical  equations  of  motion  can  now  be  obtained  by  employing 
•  Lagrange's  formulation 


ct_  _  0 

cCf 


(5) 


in  which  L  =*  Tg  +  Tp  -  V  .  Hence,  from  Eqs.  (1),  (2)  and  (3),  we  obtain 


J 

a  set  of  IT  linear,  homogeneous,  ordinary,  differential  equations'.  Such  a 
system  can  be  solved  by  assuming  harmonic  solutions  with  frequency  Cx)  of  the 
form 


%  =  \  e' 
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This,  substituted  into  equation  (6),  yields  the  set  of  linear  equations 

=°  <a) 

j 

which  has  solutions  if  and  only  if  the  determinant  of  the  coefficients  van¬ 
ishes 

^2(m-+aLj)~  kij  1=0  O) 


The  values  of  cO 6  for  which  equation  (9)  is  satisfied  are  the 

[5].  In 
here  sure 


eigenvalues  of  the  matrix  ■  I 


j  with  respect  to 


m .  .  +  a . . 
i<3  ij 


the  present  case,  in  which  both  matrices  are  positive  definite, 

N  =  3n  real  eigenvalues  and  hence  the  same  number  of  natural  frequencies* 
When  such  a  pair  of  matrices  has  been  obtained,  the  determination  of  the 
eigenvalues  is  best  performed  with  a  high-speed  computer  for  which  appropri¬ 
ate  programs  are  already  available. 

'When  the  body  is  an  elastic  continuum,  the  expressions  for  the 
kinetic  and  potential  energies  will  appear  as  integrals  and  the  system  may 
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be  considered  as  a  limiting  case  of  the  one  previously  discussed,  having  an 
inf ini. te  number  of  degrees  of  freedom.  A  procedure  for  obtaining  numerical 
solutions  for  the  eigenvalue  consists  of  approximating  the  system  by  one 
having  a  finite  number  of  degrees  of  freedom  which  can  then,  be  treated  by 
the  method  outlined  above.  This  can  be  accomplished  either  by  replacing  the 
integrals  by  quadrature  formulas  and  the  space  derivatives  of  the  displace¬ 
ments  (which  occur  in  the  expressions  for  the  potential  energy)  by  differ¬ 
ence  formulas,  or  by  variationsl  methods  such  as  those  of  Rayleigh-Ritz  and 
Galerkin. 

In  comparison  with  strip  theory,  the  foregoing  method  introduces 
the  surrounding  fluid'  into  the  problem  in  a  simple,  unified  manner,  without 
assumptions,  and  takes  into  account  quite  naturally  the  hydrodynamic  inter¬ 
ference  effects  which  are  represented  by  the  nondiagonal  elements  of  the 
jj  a^..  jj  added-mass  matrix.  In'  strip  theory  it  is  assumed  that  the  non¬ 
diagonal  elements  are  zero.  When  these  elements  are  small  in  comparison  with 
those  of  the  principal  diagonal,  application  of  the  method  of  matrix  pertur¬ 
bation  yields  a  procedure  for  correcting  the  frequencies  calculated  by  strip 
theory. 

These  considerations  sure  illustrated  in  the  following  by  several 
cases,  culminating  with  a  formulation  of  the  procedure  for  a  ship  undergoing 
a  vertical-shear  vibration. 

A  System  of  Three  Degrees  of  Freedom 

As  an  example  of  a  system  with  a  finite  number  of  degrees  of  free¬ 
dom,  consider  a  massless  string  under  tension  Z  ,  extending  from  x  -  0  to 
x  si  4 ,  to  which  small  spheres  of  radius  a  and  mass  m  are  attached  at  the 
points  x  =t  1  ,  2  ,  3  .  Let  y^  ,  y^  ,  y^  denote  small  lateral  displace¬ 
ments  in  the  y-direeticn  of  the  spheres  at  x  =  l,2,3,  respectively..  If 
this  system  is  displaced  from  equilibrium,  and  the  effect  of  the  surrounding 
fluid  is  neglected,  its  equations  of  motion  are 
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equations  (10)  "become 
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These  yield  the  secular  equation  for 
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the  solutions  of  which  are 

O ),  s*  o4  ^2  -fzT  %  <9,765  06 
6JU  as  oC  \JT  ~  I./ 1  I  4  06 
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(12) 


(13) 


Next  suppose  that  the  spheres  are  immersed  in  a  fluid  of  mass 
density  p  .  The  added  mass  of  an  isolated  sphere  of  radius  a  is  half  the 
mass  of  the  displaced  fluid.  Then,  according  to  strip  theory,  one  need  only 
replace  m  by 


M=  m+|m.p  }  mpr —irafp  (u) 

in  (10),  and  hence,  putting  fi  *  \,  Zj'iU  ,  we  again  obtain  the  expressions  (13) 
for  the  frequencies,  but  with  o £  replaced  by  &  . 
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Ibe  velocity  potential  due  to  the  motion  of  an  Isolated  sphere  with 
velocity  U  in  an  Infinite  fluid  is  that  of  a  doublet  of  strength 

/C=;-LU0?  (15) 

J  & 

at  the  center  of  the  sphere  oriented  in  the  direction  of  its  motion.  In  the 
problem  under  consideration,  however,  we  have  three  spheres  moving  in  the 
y-direction  with  velocities  y^  ,  yg  ,  and  y3  .  Since  the  radius  of  each 
sphere  is  small  in  comparison  with  the  distance  between  them,  it  is  a  good 
first  approximation  to  assume  that  the  doublet  strength  is  unaltered  by  the 
presence  of  the  other  spheres.  Ibis  approximation  is  equivalent  to  the  strip- 
theory  assumption. 

In  order  to  improve  upon  this  first  approximation  it  is  necessary  to 
take  into  account  the  image  system  within  each  sphere  of  the  two  -external 
doublets,  a  second  approximation  which  will  suffice  for  the  present  purpose. 
Tbe  image  system  of  a  y-orlented  doublet  of  strength  ^  at  a  distance  c  from 
the  sphere  is  a  similarly  oriented  doublet  at  the  Inverse  point  within  the 
sphere  of  strength yu.  a  /c  ,  and  a  line  distribution  of  doublets,  oriented  in 
the  negative  y-direction,  of  total  strength  - ^  y^a^/c3  ,  giving  a  combined 
image  strength  of  4^  a3/c3  .  Thus,  if  yg  »  1  ,  the.  doublet  strength  within 
the  second  sphere  isy^g  *  ^-a3  and  the  Bum  of  the  doublet  strengths  of  the 
image  system  ot  JU.^  within  the  first  sphere  is  J-  a6  .  According  to  the 
generalized  Taylor  added-mass  theorem  [6}  this  gives  the  induced  added  mass 


An  =irj?a6  =  |.aJrar 


In  this  manner  the  following  set  of  values  of  the  added. masses  are  obtained 


A, i  =  A2Z=  A si 

^21  —  A|2  =  Aii  =  Aa2  =■  ► 

■^31  “  A/3  =  Jg  &£  Ttfp 
Bence  the  kinetic  energy  of  the  fluid  is  given  by 

2tf - j 


(17) 


(18) 


Also  ve  have  for  the  kinetic  energy  of  the  spheres,  T  ,  and  the 

s 

potential  energy  of  the  system,  V  , 

n  T"  _ ( niZ  i  a 


2Ts=m(^  +  ^+jf] 
2V=2t(^+^  +  ^-^-^3) 


Hence,  by  (9),  the  secular  equation  for  the  frequencies  becomes 
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where  "A  * 

Since  the  terms  in  (21)  contributed  by  the  second  approximation  are 
small,  the  method  of  perturbations  is  appropriate  for  finding  the  new  eigen- 

p  p 

values.  For  this  purpose  we  first  solve  (12)  (with  oc  replaced  by  p  ) 
for  the  three  normalized  eigenvectors  ^  2,  ,  k, 3 .  ,  associated  with 


A  V/*  UiAV  MSrVli  V  V  VW*  W  — J  ^  ^  2^ 

^1  ,  (l)2  ,  CD^  .  This  gives  the  orthogonal  matrix 
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Applying  the  transformations 

we  now  obtain  from  (13),  (19),  and  (20) 
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2V  =  Z{(z-n)£^2lt 

and  the  secular  equation  "becomes 

dj\l  +  ^(|6J/I-t-l)-/31(2-v2)  0  Aol2 

o  ^(<-r)-2f  0  =0  (24) 
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Neglecting  terms  in  A  >  one  sees  that  the  eigenvalues  are  now  given  by 
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the  numerators  of  which  are  the  eigenvalues  previously  obtained  from  atrip 
theory. 

This  simple  example  illustrates  the  general  procedure  for  correcting 
the  values  of  the  frequencies  obtained  by  strip  theory  for  a  system  with  a 
finite  number  of  degrees  of  freedom. 


Next,  as  an  example  of  a  case  with  infinite  degrees  of  freedom,  con¬ 
sider  a  plate  which  bounds  an  infinite  incompressible  fluid  on  one  side.,  Tak¬ 
ing  the  x  -  plane  to  be  coincident  with  the  plate  in  its  equilibrium  posi¬ 
tion,  it  mil  be  supposed  that  the  plate  is  undergoing  a  unidirectional  oscil¬ 
lation  with  amplitude 


7^(2.  t)  -  X(x)  COS, to) t 


in  the  y-direction  at  angular  frequency  03  ,  with  nodes  at  the  points  x  *  0  , 

tr).tzir) 


*  *  * 
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The  velocity  potential  for  the  notion  of  the  fluid  in  the  region 
y  >  0  may  he  taken  of  the  form 


§  =  <p(*,y)  xi inojt 

a*  . 

(p  =  Z.  6  'iiwnx 


live  boundary  condition 


then  gives 


=  - tl  *K*A/n.riX.s--wXO‘j 


It  will  suffice  to  obtain  the  kinetic  energy  T  of  the  fluid  bounded 
by  the  planes  x  «  0  and  x  •  TT  per  unit  width  of  the  plate.  This  is  given 
by 

'fj  )j^X  ~  ji A  (*■&)  )d*d^  (28) 


where 


AM)  »  —Li  '****  nx  siin  nk  (29) 

'  U  %»9  *  9 

The  infinite  aeries  in  (29)  can  be  sunned  when  x  /  £  .  He  have  [7] 


Z.~  coz  n&  *  -Y-^t(i-2ft.6c3  0+-^,*) 


Hence,  writing  (29)  in  the  form 


we  obtain 


A(*,0  =  y  Z  ^  njz-ty-ctrs 


A(*,X)  =  'L  in  £zt<Z±ll  (3l) 
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This  shows  that  A(  x  >  %  )  has  a  logarithmic  singularity  at  x  »  ^  .  Hence, 
if  the  integral  for  the  kinetic-  energy  in  (28)  were  replaced  by  a  quadrature 
formula,  the  resulting  matrix  would  have  a  strong  principal  diagonal - 

Since  the  elements  AU  of  the  principal  diagonal  would  vary  with  i  ,  as  is 
seen  from  (31),  these  are  not  equivalent  to  a  strip- theory  added-mass  approxi¬ 
mation. 


An  exact  treatment  of  the  present  case  can  be  based  on  the  mavtm»m 
property  of  the  eigenvalues.  Put  for  the  potential  energy  per  unit  width  of 
the  plate 


r  tt 

'  (  X(ij )dadt  ■OCrSZC0t 


(32) 


and  for  its  kinetic  energy 


2Tp  =  mco* 


1  1  £(*&)  oU  ■  Amcot 


(33) 


where  (x.Jj  )  is  the  Dirac  delta  function  and  m  is  the  mass  per  unit  area 
of  the  plate.  Hence,  since  the  total  energy  is  conserved,  we  must  have 

((|w!tfAftU  +  wi.^))-k(wy}XCx)XMdjc<±^=0  (M) 

oO 

Assume  X(  x)  7L  cn  sin  nx  .  Then  (34)  becomes 
'M-t  n 


The  condition  that  cd  is  a  maximum .with  respect  to  variations  of  c1  ,  c_  . . . 

-i-  Ct 

yields 


4vr(  ^ ^ ) d/ti rn x  am n fe; ckf. cU,  =0  (36) 


a  set  of  linear  equations  for  e  ,  the  determinant  of  which  is  the  secular 

m  2 

equation  for  obtaining  the  eigenvalues  Q)  . 


"11“ 


Assuming  that  the  potential-energy  integral  is  also  diagonalized 
"by  the  transformation  X(x)  *  J~  c_  sin  nx  ,  i.e., 

vT****.  & 

>Vs;| 

frfrk(^)X(x.)XC^oLxd^  -Z^Ky, 

a  o  **> 

then  we  obtain  instead  of  (36) 


-  o 


(37) 


Hence,  setting  the  derivative  with  respect  to  c  equal  to  zero,  one  sees 
that  the  condition  for  cO  to  be  a  maximum  yields 

^  C™  +  Tn') 

or 


o)u  = 


(38) 


'Thus,  for  this  case,  the  added  mass  is  IT^  /(2n)  ,  inversely  proportional  to 
the  order  of  the  mode. 


The  Vibrating  String 

Another  illustration  of  a  case  with  infinite  degrees  of  freedom  is 
furnished  by  the  vibrating  string.  ■  It  will  be  supposed  that  the  string  of 
radius  a  is  of  infinite  length  but  restricted  to  vibration  with  nodes  at 
x  *  0. ,  ±  'TT  ,  ±2T  j  •••  along  its  length.  Let  m  denote  the  mass  of  the 
string  per  unit  length  and  V  its  tension,  and  suppose  that  the  string  is 
vibrating  transversely  in  the  y-direction  with  amplitude 

I']  (*,(:)  =■  X(x)  COS  cot  (39) 

at  angular  frequency  CO  . 

The  velocity  potential  may  be  assumed  to  be  of  the  form  [2] 

=  Cf)  0)  zinriOdt 
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where  (x,  r,  8)  are  the  cylindrical  coordinates  of  a  point,  and  (x,  r,  9) 
is  a  velocity  potential  expressible  as  the  series 


o*' 

<p  &)  *  ZI  vt-x  005  0 


(40) 


where  K^(X)  denotes  the  modified  Bessel  function  of  the  first  order.  From 
the  boundary  condition  at  the  surface  of  the  string 

(£&)  =  ir(x)an  9  t r(x)  =  -CL>X(*) 


we  obtain 


-  — ! la 

71  Vn, 


dt 


where  K^*  denotes  the  derivative  of  with  respect  to  its  argument,  and 
hence  (40)  becomes 

4  -  -  f  t  “*  (4D 

The  maximum  value  of  the  kinetic  energy  of  the  fluid  bounded  by  the  planes 
x  w  0  and  x  s  IT’  is  then  given  by 

,  ...  X  At 

2Tf  =  — pj  (^~dS  =p4)j  j  a<j>  XwdxdB 

S  '  &=0 

evaluated  on  the  surface  r  *  a  .  Substituting  from  (41),  we  obtain 

rT X  ft 


zTf  =  [  AUjH)  fe)/(^)dx  d^ 


where 


o  o 


<30 


A=  -2o?H  G(>io,)xiioinx  /W/nyd-  G(X)=-J^I 

—  AK,(/T) 


(42) 


>1  =  1 


The  modified  Bessel  function  K^{  A  )  is  defined  by  the  series  [8] 


i ,f\\_  i  i  v-  /  x  V 

K'w----i tU — ^rn^oi  (t; 


2fu 


(43) 
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Hence,  by  (42),  the  total  kinetic  energy  of  string  and  fluid  is  given  by 

2(ts+tf)  =&*“(  (48) 

and  the  condition  of  constancy  of  the  total  energy  gives 

CO1!  j  ciiij  =  rj^)2ctx  (49) 

As  in  the  previous  cn.se  of  the  vibrating  plate,  exact  solutions  for 
the  vibration  frequencies  can  be  obtained  from  the  energy  equation  (49)  by  the 
Rayleigh  -■'Hits  method.  Assume 

OO 

^  (x )  is  y  M/ro  yl  >4 

MS) 

Then  (49)  becomes  t 

cozH  <  [m-fap6rM)  = 

Ms  I  1 

The  condition  that  the  be  selected  so  that  CO  is  a  maximum  or  a  minimum 
gives 

of  r^^OLf  G  (»h)]  s;  Z  Yl 

and  hence,  denoting  the  frequency  corresponding  to  n  by  CO  ^  , 

^  I 

^1=  +  s  (so) 

where  O'  is  the  specific  gravity  of  the  string  material. 


According  to  strip  theory,  the  frequency  is  given  by 


Hence,  from  (50)  and  (51)  we  have 


<an 


(52) 
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For  harmonica  of  noderate  order  the  value  of  0(A)  in  (44)  nay  he  applied 
and  (52)  becomes 


For  example,  consider  a  steel  string  (£*■  6.8)  2  ft  long  and  0.01  ft  in 

diameter.  Since  the  foregoing  theory  is  based  on  a  string  length  of  V  be¬ 
tween  nodes,  the  value  of  a  in  (52)  will  be  taken  to  be 

a  s  &21JL  -  0.00785 

This  gives  the  following  values  of  : 


n 

G(na) 

“A. 

1 

1 

1 

10 

0.985 

1.001 

100 

0,660 

1.022 

The  example  indicates  that  the  strip- theory  approximation  would  be  an  excellent 
one  in  the  case  of  a  vibrating  string. 

Vibrating  Spheroid 

In  the  cases  of  the  Infinite  plate  and  the  Infinite  string  vibrating 
in  a  fluid,  it  has  been  possible  to  obtain  exact  solutions  for  the  natural  fre¬ 
quencies.  Since  the  principal  purpose  of  the  present  work  Is  to  illustrate  an 
approximate  numerical  procedure  for  obtaining  these  frequencies,  it  would  be 
of  interest  to  compare  the  results  of  this  numerical  procedure,  applied  to 
either  of  the  previous  cases,  with  the  known  exact  solution. 

Although  it  is  planned  to  undertake  such  a  calculation  in  the  contin¬ 
uation  of  this  study,  to  date  attention  has  been  devoted  to  the  development  of 
procedures  for  obtaining  the  added-zaass  function  for  several  cases  of  more 
practical  interest  than  the  preceding  ones  -  viz.,  a  circular  bar  of  finite 
length,  a  spheroid,  and  a  ship  form.  Of  these,  only  in  the  case  of  the  finite 
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bar,  discussed  in  Part  3  of  this  Final  Report,  has  a  complete  application 
of  the  theory,  taking  into  account  the  inertia  and  elastic  properties  of 
the  bar,  been  consummated.  The  results  obtained  in.  the  case  of  a  vibrating, 
spheroid  will  first  be  presented. 

It  has  been  shown  [9]  that  the  kinetic  energy  of  a  liquid,  for  a 
spheroid  vibrating  at  its  free  surface,  is  expressible  in  the  form 


where 

/  Ug  t  (i  iy  e  6) 


Here  , JK  ,  9)  are  the  spheroidal  coordinates  of  a  point  with  focal  dis¬ 
tance  2k  .  The  surfaces  ^  «  const,  are  prolate  spheroids,  the  value 
^  >*  £■>  q  denoting  the  equilibrium  surface  of  the  vibrating  spheroid.  The 
polar  and  equatorial  radii  of  the  spheroid  ^  are 

S  S 

The  functions  ?n  and  are  the  associated  Legendre  functions  of  the 

first  and  second  kind  [10],  and  ug  ,  vg  ,  wg  are  the  components  of  the 
velocity,  vector  vg  ,  related  to  the  velocity  vector  Ve  of  a  point  on  the 
surface  *  -l»o  *  . 

V*  =  Git  (55) 


Let  us  consider  a  pure  shear  motion  with 


^=^=0  ,  ^=4£6m.) 

Then  (54)  becomes 

a  Ttr  Tf  H  d9 


(56) 


where 


AjM 


(57) 
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Let  ua  write  the  kinetic  and  potential  energies  of  the  elastic 
spheroid  in  the  form 

zts  = 

IV  *  ]4fe^v>)X^)X(v)c^  <*>> 

where  X^O  is  the  amplitude  of  the  vibration,  given  by 

2^u.)  = -a)X(/4> 

From  the  conservation  of  the  total  energy  we  then  have 

/  Jt  Ajjs^-  kj^v)}  XjA)X(v)  cJh  ct  v>  ■  0 

where  XF 

The  Rayleigh-Ritz  nethod  will  now  be  applied,  assuming 


C*  6^ 


Then  (58)  becomes 


col  £lc.cA.-i£$h 

=  Jr,5;.  <£ 

where 

a,*  = 

^a4/  —  [( (*?  (/*,*>  6c y)dfL^ 

The  integrals  occurring  in  (59)  can  be  simplified  by  means  of  the  relations 

[8] 

fy)  —  0  '  ^ 

M«ri  =  (p„  -  pw) 
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Applylng  these  yields 


1 fy*  ?«<?■)  Pj>f> 

Bence,  sunning  with  respect  to  r,  we  obtain 

l  cA{V 

and  the  tern  containing  the  integrals  in  (59)  reduces  to 

-4X-H  (  xCy>H  (60) 

p  ■v\*»  T*h\  v.  Te-i  2e*3  / 

© 

The  condition  that  the  eigenvalues  CO  are  stationary  with  respect 
to  variations  in  the  coefficients  cr  ,  applied  to  (53)  yields  a  set  of  homo¬ 
geneous  linear  equations  and  the  secular  equation  for  determining  the  eigen¬ 
values.  The  form  of  (60),  which  represents  the  contribution  of  the  fluid  iner¬ 
tia,  indicates  that  there  will  be  several  nondiagonal  elements  due  to  the  fluid 
in  the  secular  equation.  This  will  be  studied  in  the  continuation  of  the 
present  work.  In  the  remainder  of  this  section  the  properties  of  the  added- 
mass  function  will  be  developed. 

Convergence  of  A(M.V) 

f 

For  large  values  of  n  we  have  the  asymptotic  formulas  [6] 


2ti*  l 


&V>  *-  -w  *™e) 


Hence  the  nth  term  of  the  series  for  . 


where  fl  »  cos  9  ,  V  •  cos  9*  . 

A(y*,V>)  is,  asymptotically, 

vzhr  ^owv&i)  [m  tft-ie-e')  -uk^le+e-)) 

Thus  the  series  for  converges,  except  when yU  •  V  .  Comparison  with 

(30),  the  added-mass  function  for  the  infinite  plate,  indicates  that  the  pres¬ 
ent  function  also  has  a  logarithmic  singularity  at  yK.  ■  V  and  consequently, 
if  it  were  replaced  by  a  quadrature  formula,  the  resulting  added-mass  matrix 
Ajj  would  have  a  strong  principal  diagonal.  This  Indicates  that  perturbation 
methods  should  be  suitable  for  deriving  corrected  values  of  the  frequencies 
from  those  obtained  by  atrip  theory. 


Evaluation  of 

Tables  of  values  of  and  Q^ij)  ,  which  are  required  for 

the  numerical  evaluation  of  the  added-mass  function  (57),  do  not  seem  to  be 
available  for  the  values  of  ^  near  unity  associated  with  elongated  spheroids. 
The  functions  Q^($)  may  be  successively  computed  from  the  recurrence  formu¬ 
la  (81 

n  Q^t,  (2»iti)^Qiv~(>-H)QU,  1 


where 


<&>  = 

(f-i)  4m  *  4r-0^|±i-  +2  -3^ 


Also  we  have 


of-oa;  = 


and  hence 

3L  -  (62) 

An  alternative  form,  which  will  be  useful  in  deriving  an  asymptotic  formula, 
is  obtained  by  means  of  the  relation  [8] 
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Q'ft)  = 


(?’-■)  Q,w  =  a,.,) 


These  yield  the  formula 


_  ^2.|\  ^‘»-i  (£"  0 _  (63) 

&;  ‘  £?A^Q^«-'Kfi»-,r2’>8,)+cMj’7iG),  ‘ 

From  equations  (58)  and  (59)  values  of  Q^l/Q^l  can  he  successively  computed 

It  is  of  interest  to  derive  a  series  expansion  of  Q^?)  about 
^  *  1  in  order  to  obtain  an  asymptotic  formula  for  QJVS^  for  values  of 
near  unity.  It-  will  be  shown  that 

<sy»  =  «.-h;  (e*) 


where 


a.Gr)  r 


V  >i 


LT-  -  f  ' 

rt.Ti  —  L--rr 


It  is  Interesting  to  observe  the  identity  of  the  coefficients  In 
(61)  with  those  of  the  known  expansion  [10] 


In  proving  (61)  we  will  make  use  of  the  differential  equation  satisfied  by 
the  Legendre  functions 

(•;*- ')  — -riCn-+0 

I)a§?>=0  ■ 

and  the  relation 

Q  -  Q  P-W  w/  -T-P  P 

’*0  "“ah  )  u-/~  f—  d.  fw-»  ru-A 

/ 


(66) 
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Let  us  also  write  )  as  the  Taylor  series 

W^)s=  z  y^u-o4 


*'*  *To  AJ 


Since  PQ(l)  «•  1  ,  we  immediately  have  from  (64) 


Wo)  =  H! 


From  (66)  and  the  above  differential  equation,  we  obtain 

fe20<.^2^Wn.(~7l(Yl+()Wvl.(+  =  0 

which,  by  (65)  and  (67),  can  be  expressed  as  the  series  in  terms  of  X  *  ^-1  , 

ttzz  (A-*)'.  £rrM!  £=o  rL> 

orv  p£i) 

+  2  ZI-§#  X1"  =  0 

h3o  (A-Ol 

-w 

Equating  coefficients  of  A  now  gives,  after  simplification. 


v/4 o)  -  (k-mx-™*-)  u/'t'j  pAj 

Vi  2A  M-l  A 


or,  from  (65) 


U/o )  -  M-AtOinu)  \j\i)  _  (n-A+i)'-  •  6h<> 
*’»  za>  '  t-f  /i'Z^AI 


Thus,  when  r  *»  1  ,  ve  have  by  (68) 


K(.\)  =  ^ H- 


and  similarly,  continuing  successively  with  r  *  2  ,  3  ,  ...  ,  we  obtain  ‘ 


U/%)  =  frM*0»-*frwO  uA 


It  is  now  seen  that  (61)  is  derived  by  substituting  the  expansions  (64)  and 
(67)  into  (66). 


®je  desired  asymptotic  formula  for  Q^Vo^1  is  obtained  by  substi¬ 
tuting  (64)  into  (63)  and  neglecting  powers  of  (  £  -  1)  greater  than  the 
first.  The  result  is 


-^r  |f  «  l  _l2n(n.iX0o-Hi)ti|(M)  (To) 

It  does  not  appear  to  be  useful  to  extend  this  asymptotic  formula  to 
powers  of  (^  -  1)  since,  at  the  small  values  of  n  at  which  the  formula 
gives  sufficiently  accurate  results,  better  results  are  obtained  from  the 
recurrence  relations  (61)  and  (62).  Equation  ( 70)  is  valuable  in  that  it  dis¬ 
plays  the  analytical  behavior  of  Q^/Q^  when  n  is  small  - 

At  large  values  of  n  the  iteration  procedure  for  computing  0^/$^ 
becomes  inaccurate  because  of  the  accumulation  of  errors  and  an  asymptotic 
formula  suitable  for  large  values  of  n  is  desirable.  Available  for  this  pur¬ 
pose  Is  the  expansion  [U] 


<8w  -  fa#*  ™ 

X, «  mb*, 


**Ht 


where 

v  fi  .  - 

’  }'e 
and  F  is  the  hype rgeaae trie  series 


Tbe  series  (71)  converges  when  «  1.061  ...  Wban£9<j|>^T  ,  the 

remainder  after  r  terms  of  the  series  is  less  numerically  «wn  the  (r  +  l)tb 
term;  i.e.,  the  series  is  asymptotic.  The  limiting  value  corre¬ 

sponds  to  a  spheroid  of  length-diameter  ratio  a/c  ■  3.0  ,  so  that,  for  elon¬ 
gated  spheroids,  the  asymptotic  property  of  the  series  is  very  important. 

For  m  »  1  ,  the  hypergecoetric  series  becomes 


P  -  I  _ L_  _ >1 _ I _ , 

~  Z(?n\S)  Cy-»)  ~  Q[l-nW*+S)  ‘5rrT? 


and  hence,  when  ^o<  -q-'JT ,  the  first  two  terms  of  F  will  give  a  good 
approximation  provided  the  last  term  is  very  much  less  than  unity,  Fran  the 
form  of  this  last  term  it  is  seen  that  a  nearly  equivalent  condition  is 

u+f  »prr 

But  from  the  relations 


ve  have 


cl  _  Cl 


\ _  _  a.  \ 

yrr  -  rr  “  x 

Hence  the  above  condition  becomes 


U  » 


,cu 

2  c 


-2 


(72) 


<£ 


For  sufficiently  large  n  ve  have  then 


I  + 


3  I 

_ 2ftn+0  ip 

+  2(2n+3)  |^T 


2H+2. 


jcmh^ 

3coj 


and  hence,  from  (60),  ve  obtain  after  simplification 

a«  a  afr.*-0  .  . 

(2^«)<^T+?0+ 

Values  of  calculated  by  the  exact  Iteration  formulas,  (58) 

and  (59),  and  by  the  two  asymptotic  formulas,  (70)  and  (73),  are  compared  in 
the  following  tables  for  spheroids  of  length- diameter  ratios  7.1  and  71*  In 
the  former  case  it  is  seen  that  (70)  agrees  veil  vitb  the  values  from  the  iter¬ 
ation  formulas  only  up  to  n  ■  3  ;  equation  (73)  becomes  a  good  approximation 
for  values  of  n  >  10  ,  and  indeed  yields  values  more  reliable  than  those  from 
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the  iteration  formula,  because  of  error  accumulation,  for  n  >  40  .  For  the 
much  more  elongated  body,  a/c  *  71  ,  equation  ( 70)  yields  acceptable  results 
up  to  n  *  30  .  At  n  »  70  the  value  from  (73)  is  high  by  about  4  percent, 
the  deviation  from  the  value  from  the  iteration  formula  reducing  to  about  2 
percent  at  n  «  100  . 


Q^/cy1  *  1.01  ,  a/c  as  7.14 


n 

Eqs.  (58,  59) 

Eq.  (70) 

Eq.  (73) 

1 

0.01860 

0.01857 

0.01584 

2 

0.01713 

0.01712 

0.01528 

3 

0.01572 

0.01595 

0,01446 

4 

0 » 01444 

0.01533 

0.01358 

•5 

0.01330 

0.01546 

0.01268 

6 

0.01230 

0.01649 

0.01185 

8 

0.01065 

0.01039 

10 

0.00935 

0.00920: 

12  ' 

0.00932" 

0.00822 

16 

0.00679 

" 

0.00675 

20 

'  0.00573 

0.00571  : 

30 

0.00410 

0.00409 

'  40 

0 .00322 

. ..  0  .00319 

50  . 

0.00261 

’  60 

0.00220 

70 

. 

' 

0.00191 

o 

CO 

0.00168 

* 

I 

f 

i 
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Q^1/^1  ;  ^0  *  1.0001  ,  a/c  *  71 


. 

n 

BqE.  (53,  53) 

Eq.  (70) 

Eq.  (73) 

i 

0.000199? 

0.0001997 

0.0001388 

2 

0.0001992 

0.0001992 

0.0001399 

3 

0.0001985 

0.0001985 

0.0001414 

4 

0.0001977 

0.0001977 

0.0001427 

5 

O'.  00019  69 

0.0001969 

0.0001439 

6 

0.0001958  ' 

0.0001958 

0.0001449 

8 

0.0001336 

0.0001936 

0.0001466 

10 

0.0001312 

0.0001911 

0.0001478 

12 

0.0001886 

0.0001885 

0.0001486 

16 

0.0001830 

0.0001829 

0.0001432 

20 

0.0001773 

0.0001772 

0.0001485 

30 

0.0001630 

0.0001644 

0.0001437 

40 

0.0001497 

0.0001956 

0.0001365 

50 

0.0001377 

0.0001286 

60 

0.0001272 

0.0001206 

70 

0.0001179 

0.0001131  ' 

80  ■ 

0.0001097 

0.0001070 

90 

0.0001024 

0.0000997 

100 

0.0000960 

0.0000941 

-26- 


i 

i 


i 


b 

|| 


One  final  interesting  property  of  the  added-mass  function  A(^,y) 
in  (57)  will  be  derived »  From  the  recurrence  formula  [8] 

^  +  (n+0  =  0 

we  have 


and  consequently 


£  P>  -P»>)  R,>] 

-  il^>  f?W-  -+i  fy)  (?(») 

But  we  have 

fy>=  o-^r 

P>)=  j/^0  -/T* 

and  then 

j  l#)  R  P  '^)  R  'o) 

Hence,  finally 


This  gives  the  sum  of  the  first  n  terms  of  the  series  for  A (fl,V)  w 
!-•  «l  and  n  is  chosen  sufficiently  small  so  that  Q^l/Q^l d  1  . 


(74) 
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The  equation  of  the  ship  surface  in  its  equilibrium  position  will 
be  taken  in  the  form 

,  -UTC&t 

where  x  ,  y  ,  ^  are  rectangular  Cartesian  coordinates  with  the  x  axis  in 
the  free  surface,  the  ^  axis  positive  vertically  upwards  in  the  centerplane 
of  symmetry  of  the  hull,  and  the  y  axis  completing  the  right-handed  coor¬ 
dinate  system.  From  ( 75)  one  obtains  for  the  direction  cosines  ^ 
of  the  outward  normal  to  the  surface 


1=^ 


f’+f1 

VTTS 


4=±- 


g-  -fv 


according  as  y^O  ,  where  f  *  3>f/obc  ,  etc. 

It  will  be  assumed  that  the  fluid  is  incompressible  and  inviseid, 
and  . that  the  flow  is  irrotational.  These  are  reasonable  assumptions  for  vibra¬ 
tions  of  very  small  amplitude.  Then  there  exists  a  velocity  potential  which, 
for  vibration  in  a  particular  mode  at  an  angular  frequency  <u) ,  may  be  written 
in  the  form 


The  vibration  will  be  assumed  to  be  one  of  pure  shear,  expressible  in  the  form 

v-(* ,t)  -  Vex)  eLU>t 

It  is  required,  then,  to  obtain  a  solution  of  Laplace’s  equation 

&  .  |  &<P  _  a 

3X1  +  ctW  +  U 

which  on  the  surface  (l)  satisfies  the  boundary  condition 

-fi-  =  4v«  (re) 

where  n  denotes  distance  in  the  direction  of  the  outward  normal,  and  on  the 


free  surface  the  condition 


cp  =  o 


1  BawwotaMti 
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The  amplitude  of  the  kinetic  energy  of  the  fluid  associated  with  the  vibra¬ 
tion  is  given  by 

2T=  -fj^dS  (re) 

the  integration  extending  over  the  surface  of  the  hull. 

Two  alternative  methods  for  obtaining  the  kinetic  energy  of  the 
fluid  will  now  be  formulated.  The  first  method,  in  which  it  is  assumed  that 
the  flow  is  generated  by  a  distribution  of  sources  on  the  centerplane,  requires 
the  solution  of  a  Fredholm  integral  equation  of  the  first  kind;  the  second,  in 
which  a  distribution  of  sources  on  the  hull  surface  is  sought,  leads  to  one  of 
the  second  kind.  In  principle,  an  exact  solution  for  the  latter  equation  ex¬ 
ists  and  can  be  found  by  iteration.  For  the  former,  in  general,  only  approxi¬ 
mate  solutions  can  be  found,  el  though  for  elongated  bodies  these  approxima¬ 
tions  may  be  adequate  for  practical  purposes.  Furthermore,  because  of  the 
avoidance  of  a  mathematical  singularity,  the  formulations  of  potential-flow 
problems  by  means  of  integral  equations  of  the  first,  kind  have  resulted  in 
briefer  numerical  computing  programs.  There  is  a  possibility,  however,  that 
a  procedure  for  eliminating  the  singularity  in  the  kernel  of  the  integral  equa¬ 
tion  of  the  second  kind,  which  is  presented  in  the  following,  nay  make  the  labor 
of  solving  this  type  Integral  equation  more  nearly  comparable. 

Distribution  on  Center  Plane 

The  boundary  condition  <|)  <•  0  on  the  undisturbed  position  of  the 
free  surface,  y  •*  0  ,  can  be  satisfied  at  each  Instant  by  adding  to  the  sub¬ 
merged  portion  of  the  hull  its  mirror  image  in  the  free  surface  and  treating 
the  double  form  as  a  single  body,  completely  immersed  in  an  infinite  fluid, 
and  undergoing  Vibration.  In  order  to  satisfy  the  boundary  condition  on  the 
hull  surface  (76),  assume  that  the  flow  is  generated  by  a  distribution  of 
sources  on  the  center  plane  of  the  hull  and  •  its  image,  of  strength  m^,^)  , 
where  ^  and  V  are  used  to  denote  values  of  x  and  ?  on  the  center  plane, 

4.  ” 

to  distinguish  from  the  coordinates  of  points  on  the  hull  surface. 

In  terms  of  this  distribution  the  potential  at  on  arbitrary  point 
( x  ,  y  ,  )  of  the  fluid  is  given  by 

0  =  _  J{  elyi^ 


(79) 


where 


r  =  u*-o*+f+cj-^)T 

and  the  boundary  condition  (6)  at  a  point  of  the  surface  S  yields  the 
Fredholm  integral  equation  of  the  first  kind 

»’^VW  (so) 

In  both  (79)  and  (80)  the  integration  extends  over  the  surface  of  the  center 
plane  and  its  image  in  the  free  surface.  The  Integral  equation  nay  also  be 
expressed  in  the  form 


JHvO  VoofcM) 


where 


If  (81)  were  solved,  its  solution  would  be  of  the  form 

)  =  5J *■>})$(*,})  V(*)  cU  d}  (az) 

Through  use  of  (76)  and  (79)  one  obtains  for  the  kinetic  energy  of  the  fluid 
(78) 

2TP  s  ^[jA(X',XOV(x.)V(x,<)d,xdtx.' 

where 

Al*rt  =  (a3>  - 

Here  the  hull-surface  elements  dS  occurring  in  (76)  have  been  projected  on 
the  x-£  plane,  and  only  the  projection  from  one  side  of  the  double  hull  form 
( y^  0)  has  been  taken,  since  the  energy  calculated  with  the  image  included 
is  twice  the  actual  energy.  The  function  A(x,x*)  in  (83)  is  the  desired 
quantity  from  which  on  added-mass  matrix  can  be  obtained. 


!SBSSS8B?J?.5! 


An  explicit  solution  of  (81)  can  be  obtained  for  a  very  thin  ship 
(y«  1)  .  For  this  case  the  kernel  of  (80)  behaves  like  a  Dirac  delta  function 
and  we  obtain 

*<r h.  «*> 

But  by  Cause's  flux  theorem  we  have 

<S5) 

according  as  y  £  0  ,  and,  except  near  the  borders  of  the  center  plane,  the 
corresponding  integrals  containing  derivatives  with  respect  to  y  and  |  are 
zero.  Hence,  since 

= 4  L^H* 

ws  obtain  from  (84) 

'**•(*.))  =  fjU,J)  (86) 


Comparison  of  (88)  with  the  form,  of  the  solution  assumed  In  (82)  shows  that 
hQ^,^;  x,^.)  Is  a  Dirac  S  function  which  has  the  properties  that  It  vanishes 
except  when  ^  -  x  ,  ^  ^  and  that 

f  =  _J^.  (87) 


Hence  we  obtain  for  the  added-mass  function  (83) 


I 

w~ 


ff  tfi'djdi' 
JJ 


(88) 


The  equivalent  matrix,  obtained  by  replacing  the  double  integral  by  a  quadra¬ 
ture  formula,  would  have  strong  diagonal  elements  (x  ■  x‘)  ,  although  the  non- 
diagonal  elements  would  probably  not  be  negligible. 

Better  approximations  are  obtained  by  solving  (61)  by  iteration. 
Since  k(  ^  ,  x ,  Z  )  peaks  sharply  in  the  neighborhood  of  ^  •  x  ,  ^  »  V  , 


ra^ggfHBtg  'yiUMWMMBM.1  .jJK 
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and  la  tuna'll  elsewhere,  it  is  seen  that  the  value  of  the  integral  is  affected 
only  slightly  hy  replacing  m^,^)  by  m(x,^)  in  (81);  i.e.,  the  procedure 
which  gives  an  exact  solution  when  the  kernel  is  a  Dirac  delta  function  will 
now  be  used  to  obtain  a  first  approximation.  Instead  of  (85)  we  now  employ 
the  actual  value  for  the  given  ship  fora  of  the  double  integral  of  the  kernel. 
This  yields  for  the  first  approximation 


where  , 

jl.  ^ 

and  for  the  corresponding  added-mass  function 


(S9) 


(90) 


In  order  to  obtain  a  second  approximation  we  write  (81)  in  the  fora 


-wAWJ 


=  VW  -  {J  k{t,X  iX.pd^d); 

Then,  by  the  same  reasoning  as  ws  used  to  derive  m^x ,^)  ,  we  obtain  the 
second  approximation 

=  >"•,(*, jj  -  (fj^W  -(W-11  «M|>)  <91> 

or,  from  (89), 


The  resulting  added-nass  function  is  then,  by  (62),  (85),  and  (90), 


(92) 


(93) 
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This  procedure  could,  he  continued  to  obtain  additional,  successive  approxima¬ 
tions,  but  the  second  approximation  will  probably  suffice. 


Distribution  on  Hull  Surface 


It  will  now  be  assumed  that  the  flow  is  generated  by  a  distribution 
of  sources  on  the  surface  of  the  hull  and  its  image,  of  strength  m( ^ ^ )  . 
This  gives  for  the  potential  at  an  arbitrary  point  (x,  y )  of  the  fluid 


where 


c (>  =  -[  dS 

R  =  ,  q-foss)  . 


(94) 


The  boundary  condition  (6)  at  a  point  of  the  surface  S  now  yields  the  Fred¬ 
holm  integral  equation  of  the  second  kind 


where 

k(i~) 


?  (.95) 


J 


A  procedure  for  eliminating  or  reducing  the  order  of  the  singularity 
of  the  integrand  in  (95)  when  the  points  (x,  y,^-)  and  ,  &;)  coincide 

will  now  be  indicated.  By  Gauss’s  flux  theorem  we  have 


fjlrU-'H5  =  <“> 

where  V  denotes  distance  in  the  direction  of  the  outward  normal  at  a  point 
of  the  hull  surface.  Thus  the  integral  equation  (95)  may  be  ex¬ 
pressed  in  the  form 

->n =£/<■*>  on 


which  acconrplish.es  the  desired  purpose. 
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An  approximate  solution  of  (32),  obtained  by  neglecting  the  surface 
integral,  is  given  by 


=  -£± )b*L 

>d>  HIT 

Successive  approximations  are  then  obtained  from  the  iteration  formula 

Wm;„  = 


(98) 


(99) 


which  yields  a  solution  of  the  form 


w(l'P  “  TF 


(100) 


Substituting  (76)  and  (94.)  into  (8),  we  obtain  for  the  kinetic  energy 
of  the  fluid 


zt; 


y  ( Uf  V(x> _ died*' did*' 

nJJJ  o  o 


in  .which,  by  symmetry,  only  values  of  |g  corresponding  to  y  ^  0  need  be 
taken,  the  resulting  quadruple  integral  being  multiplied  by  2  instead  of  by 
4 ,  since  the  energy  calculated  with  the  image  included  is  twice  the  actual 
energy.  Substituting  for  m(  x  ,^)  from  (100)  gives  then 

2  r  .  zeffjf n _ r  _ 

d^J  dxdx'ctj  Cty 

which  yields  the  added-mass  function 


Av,‘  Aj,,u4iiimhasMty  dj. 


(101) 


(102) 
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Stfeacagy  and  Recorsaigadations 


A  dynamical  theory  of  an  elastic  "body  vibrating  in  an  incompressible, 
inviscid  fluid  applied  in  five  examples;  system  of  three  de~ . 


grees  of  freedom,  an  infinite  plate,  &  taut  string,  a  spheroid,  and  a  ship. 

The  case  of  a  circular  bar_of  finite  length  is  separately-cons-ideretf~in~^r,t  3 
of  this  Final  Report.  *-The  first  of  these  problems  is  solved  by  employing  matr: 
theory,  the  second  and  third,  by  meaa3  of  the  method  of  Raylelgh~Ritz.  y 


^ Since  one  of  the  purposes  of  the  present  work  is  to  compare  the  Ray¬ 
leigh-hit  z  method  with  the  proposed  one  of  replacing  surface  integrals  for  the 
kinetic  and  potential  energies  by  quadrature  formulas  so  that  matrix,  methods 
become  applicable,  it  would  be  of  interest  to  obtain  solutions  for  the  same 
problem  by  both  methods .  These  can  readily  be  performed  for  either  the  infinii 
plate  or  string.  In  the  case  of  a  vibrating  spheroid,  the  mathematical  proper* 
ties  required  for  application  of  either  method  have  already  been  developed,  anc 
it  remains  to  perform  the  numerical  calculations  for  the  eigenvalues.  In  the 
last  case  a  way  of  applying  the  Rayleigh**F<its  method  is  not  apparent,  and  only 
the  formulation  for  the  matrix  method  has  been  accomplished.  „ 

Various  expressions  for  the  added-mass  function  for  a  ship,  form,  (88] 
(90),  (93),  (101),  and  (102),  have  been,  derived.  It  appears  to  be  desirable 
to  undertake  a  numerical  study  of  these  functions  for  a  particular  ship,  form, 
possibly  one  defined  by  simple  mathematical  formulas.  The  functions  defined, 
by  (88),  (90),  and  (102)  are  much  simpler  than  those  given  by  (93)  and  (101), 
although  the  latter  two  are  acre  accurate  and,  in  principle,  that  given  by 
(101)  is  exact.  A  numerical  study  would  indicate  whether  the  simpler  but  less 
accurate  forms  were  adequate  for  practical  purposes. 


would  also  be  of  essential  interest  to  compare  the  natural  fre¬ 


quencies  obtained,  with  the  added-mass  matrices  derived  from  the  above  added- 
mass  functions  with  those  computed  by  strip  theory. '  This  would  show  whether 
there  is  any  practical  advantage  over  strip  theory  of  the  present  development 
of  a  more  rational  procedure  for  the  inclusion  of  hydrodynamic  inertia  effects 
In  vibration  analyses. 
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